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We  study  analytically  the  kinetics  of  growth  of  parallel  needles,  using  a conformal  trans- 
formation to  set  up  the  iterative  nonlinear  equations.  This  allows  to  build  the  discrete 
Fokker-Planck  equation  for  the  probability  of  finding  at  time  t a given  distribution  of 
needle  lengths.  We  consider  here  two  specific  cases:  we  find  the  exact  Fokker-Planck 
equation  for  pairs  of  needles  and  its  solutions,  and  the  linear  behavior  of  a set  of  n nee- 
dles with  equal  initial  lengths.  The  corresponding  Fokker-Planck  equations  show  the 
short-wavelength  Mullins-Sekerka  instability  of  these  parallel  needles,  and  the  possible 
structure  of  the  screening  leading  to  the  scale  invariance  of  the  model. 


1 Introduction 


Laplacian  growth  is  an  ubiquitous  process  of  considerable  physical  interest:  start 
of  dendritic  growth  of  crystalline  structures  during  solidification  when  the  diffusion 
length  is  large,  electrodeposition,  viscous  fingering,  dielectric  breakdown  or  growth 
of  bacterial  colonies,  belong  to  Laplacian  growth  and  can  be  described  by  Diffusion 
Limited  Aggregation  models  (DLA)1.  Related  physical  phenomena  include  me- 
chanical cracking,  like  mud  during  drying,  crack  formations  in  pieces  of  materials 
under  strain.  The  Laplacian  field  takes  its  origin  in  the  diffusion  of  material  in 
front  of  the  growing  structures,  particles  which  aggregate,  nutrient  in  the  growth 
of  bacterial  colonies,  or  in  the  electric  field  as  in  electrodeposition  or  in  dielectric 
breakdown.  In  most  situations,  branching  occurs  during  the  growth,  leading  to 
DLA  structures.  In  spite  of  the  apparent  simplicity  of  the  physical  laws  governing 
these  systems,  and  the  very  simple  simulation  models  containing  the  essence  of  the 
physics,  our  analytic  understanding  of  DLA  remains  very  unsatisfactory.  Essence 
studies  have  used  extended  numerical  simulations  to  extract  the  scaling  laws  gov- 
erning these  systems.  Therefore  a more  modest  but  useful  approach  consists  in  a 
better  understanding  of  DLA  growth  in  the  absence  of  branching2,3.  Again,  in  this 
simpler  problem  of  growing  needles,  we  can  consider  radial  or  parallel  needles,  as 
well  as  reflection  ( model  R ) or  absorption  ( model  A4)  of  particles  (see  for  instance 
Krug5,  for  a review).  In  this  paper  we  will  concentrate  our  attention  on  systems 
of  parallel  absorbing  needles  which,  we  think,  are  of  more  fundamental  interest.  In 
two  dimensions,  conformal  mapping  allows  to  obtain  analytical  results  as  shown  by 
Derrida  and  Hakim6  for  radial  needles. 

An  important  question  in  such  systems  is  the  competition  between  the  nee- 
dles. The  instability  of  an  initial  distribution  of  equal  lengths  needles,  as  well  as 
unstable  modes  of  a growing  flat  interface  in  dendritic  growth,  is  known  as  the 
’’Mullins-Sekerka”  instability7.  In  a recent  experiment,  Losert  et  al.8  showed  the 
spatial  period-doubling  instability  of  dendritic  arrays  in  directional  solidification 
as  suggested  theoretically  by  Warren  and  Langer9.  It  will  be  the  purpose  of  this 
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preliminary  paper  to  establish  the  basic  Fokker-Planck  equation  for  the  growth  of 
n parallel  needles,  in  the  linear  regime,  which  leads  to  the  Mullins-Sekerka  behav- 
ior. We  will  indicate  in  a few  words,  on  the  basis  of  the  exact  solution  for  pairs  of 
needles,  what  we  can  expect  for  the  general  behavior  of  n parallel  needles. 

To  study  analytically  the  needle  models  of  Laplacian  growth  in  2-dimensional 
systems,  the  most  convenient  approach  is  to  use  conformal  transformations.  We  will 
only  give  in  Sec.  2,  the  main  arguments  to  set  up  these  transformations,  and  the 
reader  is  referred  to  Shraiman  and  Bensimon10,  Szep  and  Lugosi11,  Peterson  and 
Furry12,  Kurtze13,  and  Derrida  and  Hakim6.  For  the  dynamics,  we  will  follow  here 
Derrida  et  a/.6,14.  In  Sec.  3 we  write  the  closed  set  of  equations  which  describes 
the  time  evolution  in  a compact  matrix  notation.  In  Sec.  4,  we  solve  exactly  the 
general  two-needles  case,  and  establish  the  associated  Fokker-Planck  equation.  The 
case  of  n growing  parallel  needles  is  examined  in  Sec.  5,  and  the  linearized  equations 
are  derived  in  Sec.  6.  These  Fokker-Planck  equations  reveal  the  short  wavelength 
instability  leading  to  period-doubling. 

Two  ingredients  are  needed  for  the  Laplacian  growth  : on  the  one  hand  the 
Laplacian  behavior  determines  the  long  range  interaction  characteristic  of  DLA 
growth,  but  on  the  other  hand  the  inherent  Mullins-Sekerka  instability  could  not 
operate  without  the  presence  of  local  noise.  This  noise  can  be  introduced  6,14  in 
the  initial  state.  We  chose  here  to  consider  a discrete  model  with  diffusing  particles 
of  finite  size  S£  sticking  to  the  needles  at  discrete  time  intervals. 

2 The  basic  conformal  transformations 

The  classical  way  to  parametrize  in  a convenient  manner  Laplace’s  equation  with 
a zero  potential  boundary  condition  on  a set  of  n parallel  needles  is  to  introduce 
first  a transformation  which  maps  the  unit  circle  in  the  complex  plane  z onto  an 
n-branched  star  in  the  complex  plane  to: 

u>  = f{z)  = Az  n"=o  i1  ~ ei0j  /*)“’  with  E"=o  ai  = 2-  (!) 

In  this  transformation,  7roy  is  the  angle  between  two  successive  needles  {j  — 1 ,j}, 
the  sum  of  the  angles  being  27r  (equ.(l)).  The  angles  9j  fix  the  lengths  of  the 
needles,  and  A is  a parameter  which  will  be  determined  below. 

A second  transformation  maps  the  star  in  the  complex  plane  u)  into  a set  of 
parallel  needles  in  a complex  plane  : 

SI  = logw.  (2) 

In  the  plane  fl,  the  Laplacian  field  is  $ (D)  = Re  [log  (/_1  (expfl))]  . 

The  tip  positions  in  plane  z are  parametrized  by  the  angles  (pi  : 

Zi  = exp  ( icf>i ) , 0 < i < n — 1 , (3) 

and  the  lengths  of  the  needles  are  the  modulus  of  log(/(z))  at  the  points  z = Zi\ 

£i  = log  (4A)  + aj  log  |sin((0i  - 0j)/2)\  ,V  0 < i < n - 1 . (4) 

* -^j =0 
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The  n additional  constraints  must  be  imposed  onto  the  angles  fa  and  Oj  which  take 
into  account  the  fact  that  the  needle  tips  maximize  | / ( z)\  at  z = zc 

_o  otj  cot  (((/){  - Oj)/ 2)  = 0,V0<i<n  — 1.  (5) 

Now  for  the  kinetics,  the  growth  rate  is  supposed  to  be  proportional  to  the  potential 
gradient  along  the  needles.  In  addition,  the  growth  is  supposed  to  be  restricted  to 
the  tips.  Therefore,  following  refs.6,14,  where  only  the  tips  grow,  while  the  needles 
remain  at  zero  potential  ( model  A ),  the  growth  rate  of  the  needles  is: 


d£j 

dt 


oc 


=0  a3  (1  + COt2  ((&  “ 0o)/2)) 


-1/2 


,V0  <i  < n-  1 . 


(6) 


This  equation  will  be  used  to  determine  the  growth  probability  of  the  needles. 

Now  we  have  all  the  ingredients  necessary  to  build  the  recursion  relations  to 
describe  the  evolution  of  our  system. 


3 Matrix  form  of  the  basic  equations 

It  appears  first  very  convenient  to  introduce  the  matrix  C = {cij}0<i  .<n_1  where: 

dj  = cot  ((<pi  -0j)/2)  (7) 

and  condition  (5)  may  be  simply  written  (7?  is  the  vector  {cq}o<j<n-i): 

C.ct  = 0*  . (8) 

In  a similar  way,  we  can  define  D = {di0}0<i  J<n_1  where  dij  = 1 + c^.  Then 
condition  (6)  is  : 

dli/dt  oc  {D.  c?}i~1/2.  (9) 


3.1  Infinitesimal  growth 

We  consider  now  that  between  t and  t + 6t,  a particle  of  size  It  sticks  on  needle 
i.  This  increase  of  length  is  supposed  much  smaller  than  the  distance  between 
needles,  naj.  The  equations  can  be  linearized:  for  a variation  bli  of  the  length  li 
corresponds  the  variations  bfi  of  and  60 * of  Qi.  From  (4)  and  (5)  we  find: 

bli  = 6 A/ A {aj /2)  cot  {{cpi-eA/2)  60,  . (10) 

As  there  are  (2n  + 1)  variables  in  our  system  of  2n  equations  we  have  to  impose 
a supplementary  constraint.  One  choice  has  been  to  fix  one  angle  (for  example 
Qo  — 0,  ref.  14).  Here  we  will  make  the  mean  of  the  li  to  remain  zero  (this  choice 
does  not  single  out  any  angle).  This  fixes  A via  equation  (4).  Using  (10) , and 

Ei  Sii  = o. 


(6A/A)  1 = 1/(2 n)  O.C .rift  , 


(11) 


218 


where  aS&  = }0<i<Tl_1  • Here,  1 is  the  unit  matrix  and  O = {l}ij  a matrix 

filled  of  l’s.  Replacing  (11)  into  (10)  gives  6tL  as  functions  of  69j, 

SE  = -1/2  (l-(l/n)  0).C.^,  (12) 

where  (1  — (1  /n)  O)  is  a projector  on  the  subspace  orthogonal  to  the  completely 
symmetrical  component  (for  a vector  V'  it  is  Vo  = (1  /n)J2ivi)-  The  choice, 
SE[j]  = — {61  — ~}j , ...,  corresponding  to  adding  a particle  SI 

on  needle  j,  ensures  the  conservation  of  the  mean  length  (equal  to  zero).  The  t is 
for  transposition  (column  vector).  Eq.  (12)  allows  to  express  the  variations  of  69 j 
as  functions  of  the  EE.  We  can  also  relate  8<j>i  to  66 j]  Eq.  (5)  gives  : 

0 = E-Io a> {8<t>i ~ 86i) t1 + cot2 ~ • (13) 

This  makes  it  possible  now  to  calculate  the  variations  8<pi  : 

6<t>i  = (D.5^)i/(D.  at)i  . (14) 

Now,  Eq.  (9)  which  provides  the  increase  of  needle  length  per  unit  time  in  the 
diffusion  field,  enables  us  to  introduce  the  growth  probability  pi  for  a particle  of 
size  61  to  stick  on  the  tip  of  needle  i during  a time  6t.  Let  (D.W)  = If, 

Pi  = dti/CE. d ^ = (n*)_1/2  / (E  (ni)_1/2)'  (15) 

i 

The  set  of  equations,  (12) , (14)  and  (15)  together  with  (8)  define  the  successive 
iterations  of  the  growth:  if  at  time  t,  all  the  parameters  are  known,  (15)  determines 
the  growth  probabilities  on  the  tips;  knowing  the  <%  we  calculate  6 ~§  via  (12)  then 
6p  via  (14),  and  then  the  perturbation  <5C  of  C,  and  <5D  of  D.  To  calculate  6$  as 
a function  of  6Z  we  have  to  invert  Eq.  (12).  Unfortunately,  relation  (8)  shows  that 
W is  an  eigenvector  of  C for  the  eigenvalue  0,  and  therefore  C cannot  be  inverted 
directly. 

3.2  Invariants 

In  a global  rotation  86  of  the  angles  6 et  4>  the  configuration  of  the  needles  and  thus 
all  the  preceding  expressions  remain  invariant.  For  example,  C is  left  unchanged 
and  from  (12)  and  (8)  EL  also  remains  invariant: 

SE  = — (1  /n)  (1  - (1  fn)  O)  .C.  (a8$  + Stiat'j  = -(1/n)  (C  - (1/n)  O.C)  .a6&  . 

(16) 

4 The  two-needles  problem 

We  are  now  ready  to  solve  the  two  needles  case.  Here  n = 2,  but  the  needles  are 
periodically  repeated  due  to  the  cyclic  boundary  conditions  imposed  as  shown  in 
Fig.l.  With  the  zero  mean  length  constraint,  we  write, 

£o  = {Io  + h)/2  = 0;  Ji  = (4>-*i)/2. 
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Figure  1.  a)  Set  of  equally  distant  pairs  of  needles  (ao  = = 1);  b)  General  case. 


The  problem  is  completely  defined  by  the  knowledge  of  the  probability 
to  find  at  time  t a length  gap  l\.  The  general  form  of  C,  due  to  condition  (8)  and 
of  D with  elements  (dy  = 1 + cfj)  is  expressed  with  two  unknown  functions  Co(f’i) 
and  Ci(£i): 


J aic0  -a0co  \ n _ / 1 + (aic0)2  1 + (a0c0)2  \ 
\ -aiCi  Ot0c\  j ’ \ 1 + (aiCi)2  1 + (a0ci)2  / 

the  growth  probabilities  are  then  found  from  (15)  : 


(17) 


\ = 1 { y/l  + QQQlC!2]  , * 

|pi(^l)j  y/l  + a0QlC02  + yj\  + QoOlCj2  l \A  + C*oaiC02  / 

4-.1  Initial  conditions 

Fixing  arbitrarily  #o(0)  = 0,  the  initial  conditions  (i  = 0),  are: 


and  for  parameter  A, 


J <M°)  If  2arccoty/Qi/ao  ) . . 

\ 0i(O)  / \ -2arccot-v/ai/ao  J 


A = (a0)-“o/2  (ai)~ai/2  /2. 


(20) 


From  this  we  can  calculate  C (0)  and  D (0) . 
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4-2  Equations  satisfied  by  Co(£i)  and  Ci(£i) 

From  above,  we  deduce  that  the  functions  co(£i)  and  ci(£i)  verify  the  differential 
equation: 


dc0(£ i)  (X  + Qoco(^i)2)  (l  + alco(£ l)2) 

d£i  a0ai  (l  + a0aic0(£i)2)  (co(?i)  + ci(Fi)) 


and,  permuting  Co  and  c\,  a similar  equation  for  cj(£i).  The  solutions  of  these 
equations  can  be  obtained  parametrically: 

c0(Fi)  = 1 /(uv/oooi);  ci(li)  = u/yjaoai  (22) 


with  the  relation, 


The  growth  probabilities  are  then,  from  (18), 


(ao-a  i)/4 


Po(£i)  = u/(l  + u);  Pi(£i)  = 1/(1  + u). 


(23) 


(24) 


4-3  Discrete  Fokker- Planck  equation 

The  Fokker-Planck  equation  of  a pair  of  needles  is,  with  8£ i = 5£/ 2, 

P(£i,t  + St)  = p0  (Ji  - Sii'j  P{£\  - S£i,t)  +pi  (Ji  + Sii^J  P(l i + S£i,t).  (25) 

The  evolution  Eq.  (25)  can  be  expanded  to  second  order  in  6£, 

dP(l1,t)/dt  = D.dfi,P(I1,t)-v.dIi  (P{iut)  U(h))  , (26) 

where  the  following  constants  have  been  introduced, 

v = 6£/(2St)  and  D = 6£2/(8St),  (27) 

v is  the  relative  growth  velocity  of  needle  0,  D is  the  ’’diffusion”  coefficient  of  the 
sticking  between  the  two  needles,  and  £\  — ( £q  — £\)/2.  At  short  time  the  particles 
stick  at  random  on  both  needle  tips,  up  to  the  moment  when  one  needle  gives  way 
to  the  other.  The  function  U characterizes  the  screening  effect, 

W(?i)  = («(?i)-1) /(«(?i)  + 1),  (28) 

where  u{£\)  is  implicitly  defined  by  (23).  The  graph  of  U{£\)  is  shown  in  Fig. 2a. 

In  the  linear  region,  when  l\  — »■  0 , 

U(I1)^I1/(2a0a1)  (29) 

and  the  diffusion  dominates,  while  when  £\  ±oo,  U{£ i)  — » ±1,  and  the  longest 
needle  grows  at  velocity  v.  Eq.  (26)  represents  the  diffusion  of  a ’’particle”  with 
coordinate  £\  in  a potential, 


V(h)  = -v  J U(Ii) d£ i . 


(30) 
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Figure  2.  a)  Graph  of  U(t\ ) for  various  values  of  ao;  b)  Effective  potential  V ( l\ ) when  ao  = «i . 


Figure  3.  Distribution  of  probability  P{1\,  t)  at  times  t/St  = 20,30, 120. 


If  the  needles  are  equidistant,  Qo  = <*i  = 1,  and  u{i\)  = exp£1;  we  have  explicitly, 
U(t i)  = tanh  p/2).  The  potential  (shown  in  Fig.2b)  is  then  explicitly, 

(31) 


V(l\)  = — ulog  ^cosh(£i/2))  . 


The  time  evolution  of  the  solution  of  Eq.  (26)  for  equal  length  needles  (q0  = a\  = 1) 
at  t = 0 , is  shown  in  Fig.  3. 


4-4  Fluctuations  of  the  difference  between  the  branch  lengths 

The  average  of  the  difference  between  the  branch  lengths  is  zero  by  symmetry, 

(h)(t)=  f+  IiP(h,t)dI,  =0  (32) 

' ' J — oo 

because  i\P{ii,t)  is  odd.  But  the  mean  square  fluctuation 

(pi  - pi))2)  (t)  = p?)  (t)  = J_+°° IlPCh^dh 


(33) 
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is  different  from  zero.  Its  time  evolution  is  given  by  the  partial  differential  equation, 

'-dix  = (34) 


d^l)  r+°°  ^dP{iut) 

J — OO 


dt 


-oo 

/■+oo 


dt 


J °°  e\  (p.d^Pie^-v.^  (p{eut)  u{e,)))del. 

Integrating  by  parts  leads  to, 


d (l2\ 

-±^=2D  + 2v  (ixUih)). 


When  l\  <C  1,  Eq.  (35)  with  (29)  reduces  to, 

d(ll) 


dt 


2 D + 


aooti 


o2 : 

'1  / > 


(35) 


(36) 


which  has  the  following  solution  (0)  ==  O^at  small  enough  l\  or  t, 

P (~)  - l)  , (37) 

\a0aij  J 


2a0a1D 


The  growth  behaviour  is  diffusive  (^1^  ~ 2 Dt^j  , up  to  a time  tdif  — ( aoai/v ) . 

5 Mullins-Sekerka  instability  of  a comb  of  needles 

Let  us  consider  now  the  case  of  n equidistant  needles  for  which  a?  = ~ *{1, 1, 1}. 
In  this  case  it  is  convenient  to  use  the  Fourier  transform  of  the  above  equations. 


5. 1 Reduction  of  the  vector  space  dimension 

We  introduce  the  Fourier  transform  of  an  arbitrary  n-vector  V*  : 

V = F.  V with  F = {(1/n)  exp  (2?ri  jh/n)} 0<j,fc<„_i  (38) 

To  invert  Eq.  (12),  we  have  to  discard  the  completely  symmetrical  component, 
which  makes  C singular.  The  working  space  will  have  now  (n  — 1)  components,  and 
we  will  use  an  index  r to  specify  its  reduced  vectors  and  the  rectangular  matrix 
operators  relating  reduced  ( n — l)-vectors  to  n-vectors.  Thus,  we  introduce  the 
rectangular  Fourier  matrix  : 

= {-Efc/t}1<j.<n_lj0</i<n-l  • (39) 

The  corresponding  distribution  of  modes  is, 

~Lr  = Yr.~t , (40) 

a particular  mode  k,  being  given  by, 

n— 1 

£k  = (l/n)  ^ £j  exp  (27ri  jk/n) . 
j= o 


(41) 
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For  an  increase  of  length  81  on  needle  j, 

8Lr[j]  = FrM[j]  = 8I/n  {exp(27ri  ifc/n)}1<fc<7l_1  (42) 

The  reduced  Fourier  matrices  lead  to  the  following  products:  Fr.(F-1)r  = lr  and 
(F_1)r.Fr  = 1 — (l/n)0.  In  particular,  this  allows  to  write  (12)  as 

si  = — - (F_1)r.Fr.  cM.  (43) 

n 

It  is  now  easy  to  invert  Eq.  (43).  Let  Cr  = Fr.C.(F-1)r  then, 

6$  = -n  (F-1)r.(Cr)-1Fr.£2.  (44) 


5.2  Evolution  of  the  probabilities 

The  probability  P(L  ,t)  to  find  the  distribution  of  needle  lengths  L at  time  t is 
now  replaced  by  the  probability  P(Lr,  t ) to  find  a given  Fourier  mode,  and  we  have 
to  build  the  Fokker-Planck  equation  which  relates  P(Lr,  t ) and  P(Lr  + 8Lr,t-j-  8t). 

With  probability  pi  (lt  j we  add  at  time  t a particle  of  size  81  on  needle  i.  The 
corresponding  change  of  the  lengths  8t[i\  induces  a change  in  the  angles  6 and  tp, 

m = -n(F-1)ACr)-'Fr.8i\i}-,  U\i}=  {{D.^[i]}j{D.r}J-1}  (45) 


from  which  we  determine,  by  differentiation  of  Eq.  (15),  the  value  of  pi(t  + 8t): 


Pi(t  + 8t ) =pi(i)  - 


8Ui 


2H 


■3/2 


£n; 

3 


-1/2 


+ n-1/2(  £n-1/2)-2£ 


8Uh 


211 


3/2  ' 


(46) 


5.3  Initial  conditions 

In  the  present  case  for  which  the  needles  have  an  equal  length  at  t — 0,  the  initial 
angles  fa  et  8j  (0  < i,  j < n — 1)  are  also  regularly  spaced  : 

</>i(0)  = 2m /n  et  8j( 0)  = n(2j  — 1 )/n.  (47) 


5.f  Sticking  probabilities  in  the  linearized  regime 

From  expression  (45)  we  can  calculate  &J  (0)  and  <5<p(0)  corresponding  to  a variation 
brio),  for  instance  by  adding  a particle  on  needle  h.  Thus,  if  a particle  of  size  81 
is  stuck  on  needle  h of  a uniform  comb,  then  the  sticking  probability  of  the  next 
particle  is,  using  Eq.  (46),  (p(t  = 0)  = ~,Dj  = {D.T*}  , and  DJ1^2  = 1) 

• (48> 

The  sticking  probability  pi  [/),]  of  a new  particle  on  needle  i,  while  a particle  is  already 
stuck  on  needle  h ( Sih  is  the  Kronecker  symbol)  is  then, 


Pi[h]  = — + 8£ 


(n  + l)(n— 1)  csc2((i  — h)n/n) 


6 n2 


-8ih- 


2 n2 


(1  — 8ih] 


n 


(49) 
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Figure  4.  Screening  of  the  sticking  probability  by  the  longer  needle  h. 


As  expected,  in  this  result  we  see  that  the  other  needles  are  screened  by  the  longer 
needle  h.  This  is  shown  in  FigA 


6 Fokker-Planck  equation  of  the  linearized  comb 


Let  Pun(  L ,t)  be  the  linearized  probability  to  find  a set  of  needles  of  size 
L — t{£o,  i\, ...,  ih,  ■■■■>  fn-i};  the  linearization  supposes  that  the  lengths  V are 
not  too  different  in  such  a way  that  no  needle  can  be  completely  screened  (see  Sec. 
4).  In  the  linearized  approximation,  the  sticking  probabilities  on  L are  then  the 
superposition  of  the  individual  probabilities  p3  [4],  where  p3  [4]  is  deduced  from 
Pj[h]  by  replacing  8£  by  £h  (superposition  rule).  If  between  t and  t + St,  a particle 
8£  is  added  on  needle  h and  if  L'h  = I?  — si[h],  then 


n— 1 n — 1 

Plin(  L ,t  + 6t)  — EE  ^4 1 


h= 0 j= 0 


Plin{F[h],t), 


(50) 


with  V[h]  = ‘{4  + f ,h  + f ,...,4  -S£  + ^,...,4-1  + f }• 

We  recover  the  linearized  form  of  the  particular  case  n = 2 of  (25): 

Pun(£ut  + 8t)  = {p°[Q  +p0{£,1})Pun(h-Sll,t)+  [plK\  +P1[£[})Piin(£i+8£1,t) 

(51) 


as  to  first  order, 

Po  (Ji  - 8£i)  ~{1  + (£0  -h-  8V)  I A)  /2  = p°[V0}  + p 0{£[] 

Pi  (Ji  + Si i)  ~ (1  - (to  -h-  8V)/ 4)  /2  = pl\£'0]  + p^t [}.  (52) 

In  the  general  linearized  case  we  obtain  from  Eqs.  (49,50), 


PUn(£,t  + St)  = Plin(Lf[h],t)  + 

2^  ££=0  (4  - £j)  CSc2  (ih  - J>/«)  ^m(  [A],  *) 

with  L(.[/i]  = L'r—  8Lr[h]  defined  by  Eq.  (42) . 


(53) 
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6. 1 Fourier  transform  of  the  kinetic  equation,  and  the  Fokker-Planck  equation 
Equation  (53)  can  be  written  after  Fourier  transformation, 


Plin{  L 

n— 1 / 

E U- 

h=0  \ 


t + St)  = 

n— 1 

-Si  + e 

fc= i 


-2p/ifc  k 


6 n2 


fc  (1 _ i), 

n \ n ) 


P r SLr 


,t). 


(54) 


In  the  above  equation,  we  have  used  the  relation, 

^ ^ esc2  (j7r/n)  = (n2  — 6fcra  + 6fc2  — l)/3.  (55) 

Second  order  Taylor  expansion  of  the  above  Eq.  (54)  with  respect  to  SLr[h\  makes 
it  clearer.  Neglecting  the  third  order  terms  in  Si,  we  find  the  remarkable  Fokker- 
Planck  equation  (note  that  i_k  = in-k,  k being  defined  modulo  n), 

dPlin{Lr,t ) =_SiY^_d_ 
dt  6t  ^ dik 

, y^1  d2PUn(Lr,t ) Sf_Y±  d2Piin{Lr,t) 

2n26thx  dik’di -k'  2n6tkh  1 dik’dik-k- 

When  n = 2,  we  recover  (26)  linearized  with  ao  = ai  = 1,  and  U(ii)  ~ i\/2  from 
Eq.  (29). 


- (I-"  14  Plin(Lr,t ) 


d.  5 Short  discussion  of  the  linearized  kinetic  equation 
Equation  (56)  contains  three  terms  : 

i)  the  first  term  is  a drift  term,  associated  with  the  growth  velocity  v.  At  the  end 
of  the  linear  regime,  we  expect  an  exponential  screening  behavior  similar  to  U(i\) 
in  the  2-needles  case.  The  remarkable  point  lies  in  the  factor  £ (l  — ik.  If  we 
start  from  a situation  where  all  the  modes  k have  the  same  weight  ik  = i,  then  due 
to  the  ^ (l  — n ) coefficient,  the  mode  fcmax  = n/2  dominates  progressively  since  it 
has  the  highest  growth  rate.  Due  to  screening,  most  of  the  other  modes  disappear, 
and  the  system  becomes  equivalent  to  a comb  of  n/2  needles:  a succession  of  period 
doubling  is  the  result  of  the  growth. 

ii)  The  second  term  couples  the  mode  k and  n — k. 

iii)  The  third  term  is  a mode  coupling  term,  which  merges  two  modes  k'  and  k" 
into  a mode  k ~ k'  + k" . 


6. 3 Correlation  between  modes  and  fluctuations  of  a mode  q,  in  the  initial 
regime 

We  first  calculate  the  correlation  between  modes, 

W92)  (*)  = / iqiiq2P(L r,t)dL r . 


(57) 
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Integration  by  parts  using  Eq.  (56)  gives, 


9t  \£q1£q2 


lin 


- 6£  1 

(q-( 

i_^ 

\ + Sl( 

i_®^ 

91^92  & 

\ n V 

n ) 

n ) ) 

^qi^Q2  / 


(58) 


+ 


2 n28t 


kfl+92,0 


+ 


ee 

2 n5t 


(1  <5gi+Qf2,o)' 


The  fluctuation  of  a mode  q,  can  be  also  determined  from  Eq.  (56)  and  leads  to, 


dt  £ 


' 2 \ 2(51  q~ 

f ^ q 1 lin  fit  n 

Si2 


+ 


(59) 


+ 


8£ 2 


n26t 


6q,n/2  + (l  £>q,n/ 2)  1 


which  has  a solution  similar  to  the  two-needles  case,  which  underlined  the  existence 
of  a g-dependent  effective  velocity, 


veff(q)  = 


281  q 
St  n 


(-1) 


2-Di<5q|7I/2  + 2T>2  (1  — Sq,n/ 2) 


[exp  (veff(q)t)  - 1] 


' lin  veff(q) 

and  shows  that  the  mode  q = n/2  has  the  fastest  growth.  We  have  noted, 

„ 8£2  3 „ <5f2 

A = and  D2  = 


2n28t  “ 2n8t 

The  correlation  between  the  modes  has  a similar  behaviour, 

_ Dlfiqi+qifl  + -P2  (1  — ^91 +92,0) 


(60) 

(61) 

(62) 


iqiiq2/iin  veff(quq2) 

with  an  effective  velocity, 


[exp  (veff(qi,q2)t)  - 1]  (63) 


*«<«•»> -i  (2  (>-£)+ ?(*-?)) 


(64) 


Further  studies  are  in  progress  to  relate  the  linear  regime  to  the  screened  regime 
and  to  obtain  analytically  the  scaling  laws  suggested  numerically  by  Krug  et  al .4 
and  Adda  Bedia14. 
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